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ABSTRACT 


It  is  well  known  that  transient  backscatter  response  in  the  resonance  region  could  provide  useful 
information  regarding  the  natural  resonances  of  the  scattering  object,  and  such  parameters  are  often 
referred  to  as  the  singularity  expansion  method  (SEM)  parameters  (i.e.  poles  and  residues).  Poles  and 
residues,  which  respectively  represent  the  object’s  complex  natural  frequencies  and  their  corresponding 
amplitudes,  can  be  estimated  from  the  late-time  response  that  is  represented  as  a  sum  of  damped  sinusoids 
(pole  series  expansion  in  the  Laplace  domain).  Since  the  transient  backscatter  also  consists  of  the  early- 
time  response,  which  occurs  before  the  late-time,  it  cannot  be  represented  only  with  the  damped 
sinusoids.  In  the  mathematical  framework  of  the  SEM,  this  means  that  the  pole  series  converges  only  in 
the  late-time.  Here  we  consider  the  convergence  of  the  pole  series  representation  of  the  transient 
backscatter  in  terms  of  the  stability  in  the  estimated  poles  from  the  time  domain  response.  Using  a 
numerically  simulated  backscattered  response  of  a  thin  wire,  the  temporal  progression  of  the  estimated 
pole  series  is  examined  by  iteratively  estimating  them  (using  the  Matrix  Pencil  Method)  at  different  time- 
window  locations.  The  results  show  a  good  agreement  with  the  well  established  SEM  theory  (i.e.  a  stable 
set  of  poles  are  obtained  after  the  theoretically  defined  time  for  convergence).  This  iterative  approach 
could  allow  for  a  more  accurate  estimation  of  the  poles  that  represent  the  resonances  of  the  object.  We 
also  consider  the  use  of  the  “non-converging”  poles  estimated  in  the  early-time  to  approximately 
represent  the  early-time  specular  response. 


PRACTICAL  IMPLICATIONS  OF  POLE  SERIES  CONVERGENCE  AND  THE 
EARLY-TIME  IN  TRANSIENT  BACKSCATTER 


INTRODUCTION 

According  to  the  singularity  expansion  method  (SEM),  the  transient  backscattered  response  of  an 
object  can  be  represented  in  terms  of  two  components,  namely  the  early-time  and  the  late-time  responses 
[1],  The  early-time  response  is  a  forced  response  generated  while  the  incident  pulse  is  still  interacting 
with  the  object,  and  generally  contains  specular  reflections  of  the  incident  pulse  from  different  scattering 
centers  of  the  object  as  well  as  the  forced  oscillation  due  to  the  traveling  wave  on  the  object  [2],  The  late¬ 
time  response  is  generated  after  the  incident  pulse  has  passed  the  target  where  there  remains  only  the 
induced  current  oscillating  at  the  natural  resonance  frequencies  of  the  object,  and  can  be  represented  as  a 
sum  of  damping  sinusoids  (complex  exponentials),  in  which  the  complex  frequencies  correspond  to  the 
natural  resonances.  Figure  1  describes  the  transient  backscatter  response  of  an  object  in  free  space. 


Ebs(f)  1  Late-time 

v  £arly-time! - > 

«— 1 - 1 


tA=  R/c,  Time  when  the  wave  front  reaches  point  A 
tB=tA  +  Lt/c,  Time  when  the  wave  front  reaches  point  B 
tLmax=  Lmax/c,  Time  to  travel  across  the  maximum  object  dimension 
tpw=  Incident  Pulsewidth 


tL=  Minimum  time  for  convergence  of  poles  series  (late-time  resonance) 


t,  >  2t.  +  2t, 


>  2t„  +  t„ 


Figure  1.  Transient  backscatter  response  of  a  finite -size  object  in  free  space  and  the  definition  of  tL, 
which  determines  the  early-time  and  late-time  boundary.  Waveforms  launched  and  observed  from  the 

observation  point. 
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The  time-domain  SEM  representation  of  the  transient  backscattered  response  of  an  object  is 


Ebs  (it' ,  e‘ ,  es ,  t)  =  Ee  (k‘ ,  e‘ ,  es , t)u  (tL  - 1)  +  £ Rn  (k‘ , e‘ ,  es  )eA‘~h ]u  ( t  -  tL )  , 


(1) 


where  the  first  term  is  the  early-time  response  and  the  second  terms  is  the  late-time  response.  Elere  u(t)  is 
the  Eleaviside  step  function.  In  the  late -time  term,  sn=an  +  j2nfn  are  the  poles1  that  correspond  to  the 
natural  frequencies  of  the  object,  with  crn  and  fn  respectively  representing  the  damping  constants  and  the 
oscillating  frequencies,  and  Rn  are  the  residues,  which  are  the  complex  amplitudes  of  the  corresponding 
poles.  In  the  Laplace  domain,  equation  (1)  is  represented  as 


Ebs  {k‘ ,  e‘ ,  es ,  ^ )  =  Ee  (k‘ ,  g ,  es , 


^  R„{k' ,e‘  ,es)  ^V; 


s  —  s„ 


(2) 


where  r  denotes  the  Laplace  transform.  The  early-time  term  in  the  Laplace  domain  is  represented  by  an 
entire  function,  and  the  late-time  term  is  represented  by  the  expansion  of  the  poles.  The  poles  are  aspect- 
independent,  that  is,  they  are  determined  only  by  the  physical  properties  of  the  target  (shape,  size,  and 
material  property).  The  residues,  however,  are  aspect-dependent  (direction  of  incidence,  polarization  of 
incidence  and  polarization  of  scattering,  respectively  shown  as  k‘ ,  e'  and  e  )  and  related  to  the  scattering 
pattern  of  the  object.  Equations  (1)  and  (2)  corresponds  to  the  class- 1  form  of  the  SEM  representation,  in 
which  the  residue  values  are  time-independent  [3].  It  should  also  be  noted  that  a  function  continuous  at  tL, 
(rapidly  decays  afterwards)  should  better  represent  the  early-time  response.  Elowever,  the  “gated”  early- 
time  should  be  sufficient  for  the  puipose  of  describing  the  early/late-time  boundary  and  the  late-time  pole 
series  representation. 


Due  to  the  aspect-independent  poles,  which  can  be  used  as  distinctive  features  of  a  scattering  object, 
the  late-time  response  has  been  of  particular  interest  in  applications  such  as  scatterer  characterization  [4] 
and  target  identification  [5].  From  the  time-domain  backscattered  response  the  poles  can  be  extracted  by 
fitting  the  signal  as  a  sum  of  complex  exponentials  [6,7].  For  an  accurate  pole  estimation  in  practical 
applications,  a  time  window  of  the  sampled  backscattered  response  must  be  in  the  late-time  (after  tL  ,  see 
Figure  1),  where  the  signal  can  be  represented  only  by  the  second  term  in  Equation  (1).  In  the  theoretical 
framework  of  the  SEM,  this  means  that  the  pole  series  solution  is  guaranteed  to  converge  in  the  late-time 
[3].  If  the  time  window  starts  before  tL  (i.e.  the  windowed  signal  cannot  be  represented  only  by  the 
convergent  sum  of  decaying  sinusoids),  the  estimated  poles  may  also  include  the  ones  that  do  not 
contribute  to  the  convergent  pole  series,  hereafter  referred  to  as  the  non-converging  poles. 


In  this  paper,  the  pole  series  convergence  is  considered  from  a  practical  perspective  in  terms  of  the 
stability  in  the  estimated  SEM  parameters  from  the  time-domain  backscattered  response.  A  numerically 
simulated  backscattered  response  of  a  thin  wire  is  used  as  an  example.  Convergence  is  examined  by 
estimating  the  poles  at  different  locations  of  the  time  window  in  a  consecutive  manner,  which  provides  a 
temporal  progression  of  the  pole  series  representation  of  the  response.  This  approach  allows  for  more 
accurate  estimation  of  the  “true”  poles  (converging  poles  that  correspond  to  the  natural  resonances), 
without  requiring  a  priori  knowledge  of  tL,  since  they  can  be  determined  from  the  time  progression  of  the 


1  This  SEM  representation  assumes  simple,  first  order  poles  only. 
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estimated  poles.  A  similar  approach  was  used  in  the  past  to  determine  the  resonant  frequencies  of  buried 
mines  [8],  but  was  not  investigated  in  much  details.  Also  considered  in  this  paper,  is  the  use  of  the  non¬ 
converging  early  time  poles  to  represent  the  “entire  function”  corresponding  to  the  early-time  response, 
and  its  theoretical  relevance.  This  could  provide  useful  information  in  terms  of  the  impulsive  specular 
scattering  part  of  the  response,  thus  providing  possible  additional  features  for  target  characterization  and 
identification. 

SIMULATED  THIN  WIRE  BACKSCATTER 

The  transient  backscattered  response  of  a  thin  wire  model  is  simulated  using  SEMCAD,  a  finite 
difference  time  domain  computer  code  [9].  In  the  model  a  15cm  long  thin  wire  was  assumed  to  be 
perfectly  conducting  and  the  backscattered  response  was  obtained  at  2.1  meters  of  distance  by  exciting  a 
plane  wave  of  130-ps  Gaussian2  pulse.  As  indicated  in  Figure  2,  backscatter  responses  from  two  different 
wire  orientations,  i.e.  0  -  0  and  45°,  are  used.  The  computational  domain  is  truncated  with  absorbing 
boundary  condition.  The  simulated  responses  are  shown  in  Figure  3. 


Figure  2.  Thin  wire  simulation  setup.  The  incident  field  is  vertically  polarized  (0-0°),  and  the  wire  is 
parallel  to  the  plane  of  incidence.  This  configuration  allows  the  backscattered  field  to  be  co-polarized 

with  the  incident  field. 


12  14  16  18  20  22  24 

Time  (ns) 


12  14  16  18  20  22  24 

Time  (ns) 


a)  0  =  0°  b)  0  =  45° 

Figure  3.  Simulated  backscattered  thin  wire  responses 


2  It  is  practically  impossible  to  radiate  a  Gaussian  pulse  due  to  its  DC  component.  However,  it  is  used  in  the 
simulation  to  yield  a  response  that  is  close  to  the  true  impulse  response  of  the  object. 
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CONVERGENCE  OF  ESTIMATED  POLE  SERIES 

The  poles  and  the  corresponding  residues  are  estimated  using  the  matrix  pencil  method  (MPM), 
which  has  been  known  for  its  computational  efficiency  and  effectiveness  in  the  presence  of  noise[7,10]. 
In  this  process,  time -windowed  samples  of  the  response  are  represented  in  terms  of  the  sample  index  m 
and  the  sampling  interval  A,  as 

Ebs(mA()«fXzm  ,form=l,...,M,  (3) 

n= 1 

where  z  =  ev"A;  ,  M  is  the  number  of  samples  in  the  time-windowed  signal,  and  N  is  the  maximum  number 
of  desired  poles  to  estimate.  The  poles  are  estimated  by  rearranging  the  samples  into  sequential  matrices, 
and  then  solving  a  generalized  eigenvalue  problem.  Once  the  poles  are  found,  the  residues  can  be 
obtained  using  a  least  squares  approach. 

Convergence  of  the  pole  series  is  examined  by  iteratively  estimating  the  location  of  the  SEM 
parameters  in  the  Laplace  domain  at  different  time  window  locations  as  shown  in  Figure  4.  The  width  of 
the  window  (i.e.  tf-to)  was  chosen  to  be  5ns  so  as  to  include  at  least  a  few  periods  of  the  fundamental 
resonance.  The  initial  and  final  to  of  the  iterative  process  were  14.4ns  and  17.8ns,  respectively.  A  total 
of  34  different  sets  of  pole  series  were  estimated  with  the  increment  of  0.1ns.  In  Figure  5,  a  temporal 
progression  of  the  estimated  poles  (/„  only)  are  shown  as  a  function  of  to  on  the  top  of  the  original 
response  (left),  and  the  complex  poles  are  shown  in  the  Laplace  domain  (right).  In  both  cases,  it  is  shown 
that  the  poles  converge  to  a  stable  set  of  locations  after  tL  ,  which  shows  a  good  agreement  with  the  SEM 
framework,  where  tL  is  defined  as  the  minimum  time  for  the  pole  series  convergence  [3].  It  is  also 
observed  that  only  the  symmetric  mode  resonances  are  extracted  for  0  =  0° ,  whereas  both  symmetric  and 
non-symmetric  modes  are  extracted  for  0  =  45“ .  The  non-converging  poles  that  are  extracted  before  tL 
vary  in  time  and  do  not  represent  the  natural  resonances  of  the  scatterer,  since  their  values  are  the  results 
of  attempt  to  fit  the  early-time  signal  to  the  late-time  representation.  Some  of  these  poles  may  be  the 
result  of  the  MPM’s  attempt  to  estimate  the  poles  in  the  second  layer  [1],  which  are  very  difficult  to  excite 
and  estimate.  Figure  6  also  shows  the  temporal  progression  of  the  poles  in  terms  of  their  corresponding 
residue  amplitudes.  In  the  figure,  it  is  clearly  seen  that  the  residues  of  the  non-converging  poles  vary  in 
time  more  than  that  of  the  converging  poles  while  that  of  the  converging  poles  are  relatively  constant. 


Figure  4.  Iterative  pole  extraction  approach  by  sliding  the  sample  time  window 


Amplitude  Amplitude 
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(b)  9  =  45° 

Figure  5.  Estimated  poles  from  the  wire  backscatter  responses.  The  plots  on  the  left  show  the  poles 
(frequencies  only,  no  damping  constants)  estimated  at  corresponding  time  window  starting  points  as  the 
window  is  shifted.  The  plots  on  the  right  show  the  same  poles  in  the  Laplace  domain  (damping  constant 
vs.  frequency).  Convergent  poles  show  stability  in  their  values  with  time  window  shift. 


Time  (ns) 


Time  (ns) 


Figure  6.  Estimated  poles  (frequencies  only,  no  damping  constants)  plotted  in  terms  of  corresponding 


residue  amplitudes,  \Rn  .  (log  scale) 
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Since  the  converging  poles  can  be  clearly  determined  from  the  temporal  progression  in  Figure  5,  this 
approach  could  be  an  effective  way  to  accurately  estimate  resonant  parameters  without  requiring  an  initial 
guess  for  the  early-time/late-time  boundary,  especially  if  the  scatterer  is  unknown. 


REPRESENTATION  OF  THE  EARLY-TIME  USING  ESTIMATED  POLE  SERIES 

Another  useful  insight  is  made  regarding  the  non-converging  poles  when  ta  is  in  the  early-time 
(before  tL  ).  Some  additional  target  characteristics  may  be  examining  the  non-converging  poles.  Figure  7 
shows  the  estimated  poles  for  0  =  0°  when  t0  =  14.74ns.  From  the  temporal  progression  plot  in  Figure  5, 
it  is  possible  to  distinguish  the  converging  poles  from  non-converging  ones  based  on  the  consistency  of 
their  locations  in  the  Laplace  domain  as  shown  in  Figure  7.  Therefore  this  pole  series  can  be  separated 
into  two  different  sums,  namely  np  series  and  p  series,  corresponding  to  the  non-converging  and 
converging  poles,  respectively  as  follows: 

E*  (<;>„  -  14.74ns)  =  2  -  O  +  Z  "  O  •  (4) 

« P  P 

Using  equation  (4),  the  response  is  reconstructed  into  two  different  signals  as  shown  in  Figure  8.  The 
reconstructed  signal  from  p  series  corresponds  to  the  slowly-varying  part  in  the  early-time  and  continues 
into  the  late-time,  closely  matching  the  original  response.  The  signal  reconstructed  from  np  series  (dashed 
green  line)  is  more  impulsive  than  that  reconstructed  from  p  series  (solid  blue  line)  and  vanishes  after  tL, 
which  indicates  that  np  series  contributes  to  the  early-time  forced  response. 
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Figure  7.  Estimated  Poles  for  0  =  0°  when  (  =1 4.74ns.  Using  the  data  from  Figure  5a,  the  converging  and 

non-converging  poles  can  be  determined. 


Practical  Implications  of  Pole  Series  Convergence  and  the  Early-time  in  Transient  Backscatter 


7 


Figure  8.  Reconstructed  signals  from  two  separate  pole  series  (estimated  using  t0=14.74ns):  converging 
poles  (solid  blue),  non-converging  poles  (dashed  green).  Original  response  is  shown  in  dotted  red. 


In  the  Laplace  domain  equation  (4)  can  be  expressed  as3 


R 


np 


nP  $  Snp 


-I 


R„ 


p  s~sp 


(5) 


According  to  the  SEM  framework,  the  early-time  representation  should  be  an  entire  function,  which  means 
that  the  first  term  in  equation  (5)  is  not  a  valid  representation.  However,  we  can  relate  the  term  to  the 
generalized  representation  of  an  entire  function  in  the  left-half  plane  using  the  Cauchy  integral  [11],  which 
is 
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ds'  ,  where  c 


lim  c  (left-plane  contour). 
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(6) 


The  Cauchy  integral  in  equation  (6)  can  be  considered  as  poles  with  continuous  residue  distribution 
R(s')  at  s  — >  oo  ,  suggesting  that  the  entire  function  can  be  regarded  as  a  function  with  the  singularities  at 
infinity.  Therefore,  from  this  perspective  the  pole  (or  singularities)  in  the  first  term  of  equation  (5)  can  be 
approximately  represented  as  an  “entire  function”  to  describe  the  early-time  forced  response,  which  could 
provide  useful  information  in  the  impulsive  part  of  the  response,  thus  providing  possible  additional 
features  for  target  characterization  and  identification. 


Shifted  left  in  time  by  t0  for  the  simplicity  of  the  expression  without  loss  of  generality 
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SUMMARY  AND  CONCLUSION 

In  this  paper,  the  convergence  of  the  SEM  pole  series  in  transient  backscatter  was  explored  in  terms 
of  the  stability  in  the  estimated  SEM  parameters.  Using  the  simulated  response  of  a  thin  wire,  the  pole 
series  convergence  was  investigated  through  an  iterative  approach  that  provided  the  temporal  progression 
of  the  estimated  poles.  The  analysis  showed  good  agreement  with  the  SEM  theory,  i.e.  the  pole  series 
convergence  was  observed  after  the  theoretically  defined  tL ,  the  minimum  time  for  convergence.  The 
iterative  pole  extraction  approach  was  also  shown  to  be  an  effective  method  of  accurately  estimating  SEM 
parameters,  particularly  when  the  scatterer  is  unknown. 

It  was  also  suggested  that  the  use  of  the  non-converging  early-time  poles  could  approximately 
represent  an  “entire  function”  to  describe  the  early-time  forced  response.  This  response  could  provide 
additional  useful  information  about  an  unknown  scattering  object,  e.g.  the  features  extracted  from  the 
early-time  forced  response  could  be  used  to  better  characterize  or  identify  unknown  scatterers. 

This  report  only  considered  a  simple  object  in  a  noiseless  environment  (simulated).  It  would  be  of  an 
interest  to  investigate  various  types  of  objects  in  more  complex  environment  (e.g.  noisy  responses, 
multiple  scatters,  etc.)  for  more  practical  examples. 


Practical  Implications  of  Pole  Series  Convergence  and  the  Early-time  in  Transient  Backscatter 
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